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in  the  technical  literature.  A second  report  (Reference  1-2)  docu- 
ments this  effort.  The  emphasis  was  on  the  upper  stratosphere,  meso-. 
sphere,  and  lower  thermosphere,  where  most  HF  wave  absorption  occurs. 

The  minor  species  included 

(0) , atomic  oxygen 
(0_) , ozone 

(02(a  Ag)),  molecular  oxygen,  excited  a-state 

(N) , atomic  nitrogen 

(NO) , nitric  oxide 

(N02) , nitrogen  dioxide 

(H20) , water  vapor 

(C02) , carbon  dioxide. 

Many  other  species  are  important  to  the  general  atmospheric  model,  but 
these  are  expected  to  have  the  primary  direct  effect  on  the  ionization 
level,  and  therefore  on  the  electromagnetic  wave  absorption. 

METHOD 

The  present  report  documents  our  findings  on  the  sensitivity  of  the 

HF  wave  absorption  to  the  expected  variations  in  these  minor  atmospheric 

species.  A summary  of  conclusions  is  listed  below,  but  first  we  pre- 

% 

sent  a brief  description  of  the  tools  employed  in  the  analysis  and  the 
organization  of  the  rest  of  the  report. 

The  basic  analysis  tool  is  the  WEPH  VI/ROSCOE  D-region  atmospheric 
chemistry  model  (References  1-3  and  1-4).  The  model  is  a considerable 


Reference  1-2.  Jordano,  R.J.,  et  al,  Suggested  Natural  Variations  in 
Atmospheric  Minor  Neutral  Species , GE78TMP-26  (ARBRL-CR-00375) , 
General  Electric  Company — TEMPO,  July  1978. 

Reference  1-3.  Knapp,  W.S.,  et  al,  WEPH  VI:  A Fortran  Code  for  the 
Calculation  of  Ionization  and  Electromagnetic  Propagation  Effects 
Due  to  Nuclear  Detonations  ( U)t  Volume  3:  Computational  Models  (V), 
GE75TMP-53  (DNA  3766T-3) , General  Electric  Company— TEMPO,  October 
1976  (Confidcntial/Formcrly  Restricted  Data). 

Reference  1-4.  Knapp,  W.S.,  et  al.  The  ROSCOE  Manual:  Volume  II:  At- 
mospheric Chemistry  Models , GE74TMP-59  (DNA  3964F-11),  General 
Electric  Company — TEMPO  (to  be  published). 
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simplification  from  the  typical  multi-species  chemistry  models  that 
solve  the  coupled  non-linear  differential  rate  equations  by  numerical 
integration.  The  portion  of  the  WEPII/ROSCOE  model  that  we  use  is  the 
steady-state  ion  species  model  which  has  shown  particularly  good  agree- 
ment with  the  more  detailed  models  in  the  altitude  regime  where  HF  ab- 
sorption takes  place  (between  SS  and  115  km  altitude).  The  WEPH/ROSCOE 
model  is  of  a form  that  is  conducive  to  repetitive  parametric  sensi- 
tivity analyses.  Hopefully,  important  conclusions  can  be  verified  with 
the  more  detailed  (and  more  expensive)  numerical  models.  The  WEP11/ROSCOE 
model  accepts  the  neutral  species  concentrations,  and  the  various  ion 
production  rates,  as  inputs.  It  then  computes  steady-state  concentra- 
tions of  electrons  and  a few  dozen  ion  species.  For  the  levels  of  nu- 
clear disturbance  investigated  in  this  analysis,  the  energy  deposition 
rates  arc  sufficiently  small  that  the  neutral  species  concentrations 
are  not  disturbed.  The  species,  reactions,  and  method  of  solution  of 
the  model  are  summarized  in  Appendix  A. 

A second  analysis  tool  employed  in  the  present  study  is  the  natural 
ion  production  rate  model.  The  solar  flux  impinging  on  top  of  the  at- 
mosphere and  the  interaction  cross-sections  with  the  air  species  are 
complicated  functions  of  photon  wavelength.  Computation  of  the  solar 
ion  production  rate  at  any  point  in  the  atmosphere  requires  numerical 
integrations  over  these  spectral  functions.  We  have  assembled  simple 
analytic  representations  of  the  natural  solar  and  cosmic  ion  production 
rates  in  which  much  of  the  spectral  integration  is  precomputed.  The 
formulations  do  retain  the  necessary  functional  relationships  to  reflect 
proper  sensitivity  to  variations  in  the  atmospheric  parameters,  however. 
These  model  formulations  are  presented  in  Appendix  B. 

Section  2 of  this  report  presents  a general  analysis  of  atmospheric 
IIF  wave  absorption,  both  in  the  natural  ionosphere  and  under  conditions 
of  nuclear  disturbance.  This  is  intended  as  a general  orientation,  in 
preparation  for  the  sensitivity  analysis  which  is  presented  in  Section 
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SUMMARY  OF  CONCLUSIONS 

This  section  presents  a brief  summary  of  conclusions  obtained  from 
the  sensitivity  analysis 

• Daytime  high-frequency  radiowave  absorption  occurs 
between  55  and  115  km  in  the  natural  ionosphere. 

Typically  two-thirds  of  the  absorption  is  caused  by 
the  Lyman-alpha  radiation  in  the  mesosphere,  with 
the  remainder  distributed  between  cosmic  ray  (upper 
stratosphere)  and  x-ray  and  ultraviolet  (lower  thermo- 
sphere) sources . 

• During  the  recovery  phase  of  an  HF  circuit  following 
nuclear  blackout,  the  nuclear  enhanced  absorption  due 
to  delayed  radiations  occurs  below  80  km  altitude. 

• In  terms  of  a lumped-parameter  solution,  the  electron 
concentration  may  be  represented  by  the  following  ex- 
pression in  the  regions  where  significant  HF  wave 
absorption  occurs: 

(e)  21  dTa  N^d  cra’3 
where 

-3  -1 

q = total  ion  production  rate  (cm  s ) 

A = total  electron  attachment  rate  (s  *j 

D = effective  electron  detachment  rate,  depending 
on  distribution  of  negative  ions  (s  *) 

a,  = effective  electron  recombination  rate,  depending 

d -3  -1 

on  distribution  of  positive  ions  (cm  s ). 

• Increasing  the  minor  oxygen  species  concentrations  (0), 

(03) , (02(a1Ag)),  generally  increases  the  integrated 
vertical  one-way  HF  absorption  through  the  atmosphere. 
Species  variations  of  a factor  of  9 give  about  a factor 
of  1.25  variation  in  absorption  in  the  natural  ionosphere, 
and  typically  a factor  of  2 variation  in  the  nuclear  dis- 
turbed ionosphere. 
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• Variation  of  the  oxygen  species  has  impact  generally 
through  the  detachment  rate  D.  Both  (0)  and  ((^(a^Ag)) 
enhance  the  detachment  rates  for  0^  and  0".  Increasing 
(0)  also  tends  to  increase  the  relative  proportions  of 
(02)  and  (O'),  while  increasing  (0^)  and  (02(a1Ag))  has 
the  reverse  effect.  Thus  (0)*is  the  most  important 
minor  oxygen  species. 

• Increasing  the  minor  nitrogen  species  concentrations  (N) , 
(NO),  (N02)  generally  increases  the  integrated  vertical 
one-way  HF  wave  absorption  through  the  atmosphere.  Species 
variation  of  a factor  of  9 give  about  a factor  of  2 varia- 
tion in  absorption  in  the  natural  ionosphere,  and  typically 
less  than  a factor  of  1.1  variation  in  the  nuclear  dis- 
turbed ionosphere. 

• The  most  important  nitrogen  species  is  (NO)  which  acts 
primarily  through  the  (Lyman -alpha)  ion  production  rate  q. 

• Increasing  (CO.,)  and  (H20)  generally  decreases  the  inte- 
grated vertical  one-way  HF  wave  absorption  through  the 
atmosphere.  Species  variations  of  a factor  of  9 give  only 
a factor  of  1.1  to  1.2  variation  in  the  absorption,  acting 
mainly  through  . 

RECOMMENDATIONS 

The  important  species  identified  are  (0),  (0^) , ((^(a^Ag)),  and 
(NO),  with  the  oxygen  species  being  most  important  under  nuclear  dis- 
turbed conditions  and  the  nitric  oxide  being  most  important  under 
natural  conditions.  Variations  of  an  order  of  magnitudfc  in  the  daytime 
concentrations  of  these  species  are  certainly  possible  as  a function 
of  latitude  and  season.  This  could  yield  a factor  of  2 variation  in 
the  integrated  one-way  vertical  HF  wave  absorption.  This  is  comparable 
to  the  potential  variations  due  to  changes  in  temperature,  pressure, 
and  mass  density  that  were  observed  in  an  eralier  study.  These  varia- 
tions arc  important  to  the  prediction  of  HF  circuit  recovery  following 
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nuclear  blackout  and  the  variable  atmospheres  should  be  incorporated 
into  the  nuclear  effects  computer  models.  Additional  atmospheric 
measurements  are  required  to  expand  and  verify  our  understanding  of 
the  minor  species  variations,  particularly  at  mesospheric  altitudes. 
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SEU1UN  2 

BASELINE  ANALYSIS 
ATMOSPHERIC  HF  WAVE  ABSORPTION 


NATURAL  IONOSPHERE 

The  total  natural  ion  production  rate  for  mean  conditions  at  noon 
(sunspot  number  = 60,  latitude  = 45  decrees,  spring/fall  seson,  solar 
zenith  = 45  degrees)  is  shown  in  Figure  2-1.  Note  that  the  cosmic  ray 
source  dominates  from  the  surface  up  to  about  65  km  altitude,  the  Lyman- 
alpha  source  dominates  between  65  and  95  km  altitude,  and  the  UV  source 
dominates  above  95  km  altitude.  The  solar  x-ray  source  makes  a second- 
ary contribution  above  about  90  km.  The  production  rate  profile,  of 
course,  changes  for  other  conditions,  but  the  essential  character  of 
the  profile  remains  the  same  as  in  Figure  2-1.  We  shall  discuss  each 
of  the  kinds  of  variations  in  the  following  paragraphs. 

The  variations  with  solar  activity  are  easy  to  visualize  from  Fig- 
ure 2-1.  The  cosmic  ray  source  has  a total  variation  of  about  a factor 
of  3:  that  is,  an  increase  of  about  a factor  1.5  over  that  in  Figure 
2-1  for  solar  minimum  conditions,  or  a decrease  of  a factor  of  0.5  from 
that  in  Figure  2-1  for  solar  maximum  conditions.  Note  the  inverse 
dependence.  The  Lyman-alpha  and  UV  sources  show  a total  variation  of 
about  a factor  of  2 (decreasing  a factor  of  0.67  or  increasing  a factor 
of  1.33  from  that  in  Figure  2-1,  for  solar  minimum  or  maximum  conditions, 
respectively).  The  x-ray  source  varies  the  most.  It  may  decrease  a 
factor  of  10  or  more  from  that  in  Figure  2-1  for  solar  minimum  condi- 
tions, or  increase  a factor  of  2 for  solar  maximum  conditions. 

The  variations  with  latitude,  season,  and  time  of  day.  all  have  a 
strong  variation  characterized  by  the  solar  zenith  angle,  plus  some 
secondary  variations  due  to  changes  in  the  atmospheric  pressure,  density, 
and  major  and  minor  species  concentration  profiles.  We  shall  discuss 
the  variations  with  solar  zenith  first.  The  cosmic  ray  source  is 
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Figure  2-1.  Mean  daytime  natural  ion  production  rate  profile 
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independent  of  solar  zenith.  Each  of  the  other  sources  has  the  zenith 
angle  dependence  in  the  attenuation  function  which  determines  how  deep 
into  the  atmosphere  the  radiation  can  penetrate.  Thus,  the  variation 
in  solar  zenith  effectively  raises  or  lowers  the  underside  of  the  given 
source  profile.  From  the  secy  dependence,  one  can  show  that  the  alti- 
tude for  a given  number  of  attenuation  e-foldings,  or  the  effective 
altitude  of  the  profile  underside,  varies  as 

h = h0  " Hs  ^n(cosx)  (2-1) 

where  hQ  is  the  altitude  for  an  overhead  sun  (x=0),  and  Hs  = pressure 
scale  height  (e-folding  distance).  The  pressure  scale  height  in  the 
underside  regions  (60-85  km  for  Lyman-alpha,  80-100  km  for  x-ray,  and 
90-100  km  for  UV  sources)  is  about  6 to  7 km.  Thus  for  example,  the 
effective  height  of  each  of  these  profiles  in  Figure  2-1  may  be  lowered 
about  2 km  for  an  overhead  sun,  or  may  be  raised  about  2 km  for  x = 60 
degrees,  raised  about  6 km  for  x = 7$  degrees,  or  raised  about  9 km 
for  x = 80  degrees.  (For  higher  zenith  angles  the  secy  dependence 
breaks  down.) 

In  addition  to  the  solar  zenith  angle  dependence,  there  are  some 
secondary  variations  in  the  ion  production  rate  due  to  atmospheric 
variations.  The  cosmic  ray  source  is  proportional  to  the  total  number 
density  which  can  have  seasonal  variations  of  10  percent  at  mid-lati- 
tudes, and  40  percent  at  the  higher  latitudes  (the  densities  are  low 
in  winter  and  high  in  summer  in  the  important  altitude  regime).  The 
Lyman-alpha  source  is  proportional  to  the  nitric  oxide  concentrations. 
Seasonal  variations  of  (NO)  of  plus  or  minus  a factor  of  2 (maximum  in 
summer)  are  expected  at  mid- latitudes,  and  higher  concentrations  (factor 
of  4)  at  all  seasons  arc  expected  at  higher  latitudes.  The  attenuation 
of  the  Lyman-alpha  source  is  determined  by  the  line  integral  of  (02) 
penetrated,  which  varies  approximately  like  the  local  pressure.  At 
the  altitudes  of  interest  to  Lyman-alpha  attenuation,  the  pressure 
varies  about  10  percent  (maximum  in  summer)  at  mid-latitudes  and  about 
40  percent  at  higher  latitudes.  At  the  higher  latitudes,  this  effectively 
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raises  (in  sununcr)  or  lowers  (in  winter)  the  altitude  of  the  undersiu. 


of  the  Lyman-alpha  profile  by  a few  kilometers.  Similar  variations 
will  be  seen  in  the  x-ray  and  UV  source  profiles.  In  the  90  to  100  km 
region,  these  profiles  will  be  effectively  raised  or  lowered  slightly, 
in  the  summer  or  winter  seasons,  respectively,  due  to  variations  in 
the  major  species  concentrations.  The  UV  ionization  of  atomic  oxygen 
will  be  significantly  reduced  at  higher  latitudes,  due  to  the  lower 
levels  (factor  of  10)  of  (0).  A couple  of  examples  of  natural  varia- 
tions in  the  daytime  ion  production  rate  are  shown  in  Figure  2-2. 

The  total  ion  production  rate  of  Figure  2-1  is  partitioned  into 
particular  ion  production  rates  according  to  the  prescriptions  defined 


in  Equations  B-4,  B-24  to  B-26,  and  B-31  through  B-35.  These  phase  one 
ion  production  rates  are  shown  in  Figure  2-3.  The  N^»  N*,  and  0*  under- 
go rapid  charge  transfer  to  0*  and  NO*,  as  described  in  Equations  A-l 
and  A-2,  and  the  effective  0*  and  N0+, phase  two  ion  production  rate 
profiles  used  in  the  chemistry  model  are  shown  in  Figure  2-4.  Note 
that  the  dominant  ion  production  rate  is  for  NO*  between  about  65  and 
95  km,  due  primarily  to  the  Lyman-alpha  source.  Below  65  km,  the  domi- 
nant production  rate  is  for  0*,  due  primarily  to  charge  transfer  from 
the  N^  produced  by  cosmic  rays.  Above  about  9S  km,  the  dominant  pro* 
duction  rate  is  also  for  02>  due  primarily  to  direct  UV  ionization. 

The  resulting  steady-state  total  ion  and  electron  concentrations 
for  this  daytime  case  are  shown  in  Figure  2-5.  The  HF  wave  absorption 
is  due  primarily  to  electron-neutral  collisions,  and  is  computed  fzum 


where 


4.6  x 10  (e)v 
u)2  + (0.775V)2 


dB/km 


(2-2) 


(e)  = electron  concentration  (cm*  ) 
u>  = wave  carrier  frequency  (radian  s'1) 
v = electron  neutral  collision  frequency  (s*  ). 


The  collision  frequency  v is  computed  from 
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Figure  2-2.  Alternate  daytime  natural  ion  production  rate  profil 


le  natural  phase  1 Ion  production  rate  profiles. 


Figure  2-4.  Mean  daytime  natural  phase  2 ion  production  rate  profiles. 


SPRING/FALL 
LATITUDE  = 
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v=  2.6*  10"n(N2)T  + l.S  x 10_11(02)T+  8.2  * l(f10(0)T*i  s*1 

(2-3) 

where  (N^) , (C^) , and  (0)  are  species  concentrations  (cm-3),  and  T is 
the  electron  temperature  (K) . Note  that  below  about  90  km  altitude, 
the  third  term  may  be  neglected  and  the  expression  for  v simplifies 
to 

v = 1.7  X 105P  s"1  (2-4) 

_2 

where  P = atmospheric  pressure  (dyne  cm  ).  The  electron  neutral  col- 
ision  frequency  is  also  shown  in  Figure  2-5.  Note  that  the  actual 
quantity  plotted  is  v/(2ttx10^)  which  gives  the  collision  frequency 
in  units  directly  comparable  to  the  carrier  frequency  in  megahertz. 

Thus  for  an  HF  carrier  frequency  of  10  MHz,  for  example,  w is  less  than 
V below  about  55  km  altitude  and  the  absorption  is  approximately  propor- 
tional to  (e)/v.  Above  about  55  km  altitude,  w is  greater  than  v and 

2 

the  absorption  is  approximately  proportional  to  (e)v/w  . For  the  typi- 
cal HF  band,  this  crossover  altitude  is  seen  to  vary  from  about  65  km 
(for  3 MHz)  to  about  45  km  (30  MHz) . 

The  incremental  and  cumulative  (one-way  vertical)  absorption  for 
a 10  MHz  signal  are  shown  in  Figure  2-6.  Note  that  99  percent  of  the 
absorption  occurs  between  55  and  115  km  altitude.  The  three  bumps  in 
the  incremental  absorption  profile  correspond  to  regions  where  different 
ionization  sources  dominate.  The  bulk  of  the  absorption  (about  60  per- 
cent) occurs  in  the  central  bump  between  about  70  and  95  km  altitude, 
where  the  Lyman-alpha  source  dominates.  About  25  percent  of  the  ab- 
sorption occurs  in  the  upper  bump  above  95  km,  where  the  UV  and  x-ray 
sources  dominate.  Only  about  15  percent  of  the  absorption  occurs  below 
65  km,  where  the  cosmic  ray  source  dominates. 

These  observations  change  somewhat  for  other  conditions.  Figure 
2-7  compares  incremental  absorption  profiles  for  the  alternate  ioniza- 
tion profiles  of  Figure  2-2.  In  all  cases,  the  absorption  occurs 
between  55  and  115  km  altitude.  For  the  midlatitude  summer  profile. 


INCREMENTAL 

ABSORPTION 


INCREMENTAL  ABSORPTION  (dB/km)  OR  CUMULATIVE  ABSORPTION  (dB) 
Figure  2-6.  Mean  daytime  natural  wave  absorption  at  10  MHz. 


the  UV,  x-ray,  and  Lyman-alpha  energy  all  penetrate  deeper  into  the 
atmosphere,  tending  to  swamp  the  contribution  of  cosmic  rays  to  the 
absorption.  The  proportional  contributions  to  the  total  absorption 
become:  UV  and  x-rays:  34  percent;  Lyman-alpha:  66  percent;  cosmic  rays: 
less  than  1 percent.  For  the  high-latitude  winter  profile,  the  larger 
atmospheric  attenuation  of  the  solar  flux  severely  reduces  the  UV  and 
x-ray  contributions,  and  effectively  raises  the  altitude  band  of  Lyman- 
alpha  influence  to  between  75  and  105  km.  The  cosmic  rays  become  rela- 
tively more  important  and  the  proportional  contributions  to  the  total 
absorption  become:  UV  and  x-rays:  3 percent;  Lyman-alpha:  65  percent; 
cosmic  rays:  32  percent.  Thus,  it  appears  that  for  most  daytime  cases 
of  interest,  the  Lyman-alpha  ionization  contributes  about  two-thirds 
of  the  total  HF  absorption  and  any  part  of  the  remaining  one-third  may 
be  contributed  by  cosmic  rays,  or  solar  x-rays  and  UV. 

The  effective  lumped-parameter  reaction  rate  coefficients  are  shown 
in  Figure  2-8.  The  attachment  and  detachment  models  are  described  in 
Appendix  A.  Below  90  km  (and  throughout  the  region  where  attachment 
is  important) , the  dominant  attachment  processes  are  the  three  body 
reactions  (Reactions  #166  and  168)  that  involve  and  C^,  and  form  C^- 
Above  45  km  (and  throughout  the  region  where  detachment  is  important) , 
the  dominant  detachment  processes  are  the  collisional  detachment  (Reac- 
tion #177)  and  associative  detachment  (Reactions  #173,  174,  and  176). 
These  reactions  involve  (^(^A)  and  0,  and  remove  electrons  from  0^  and 
O-.  From  Equation  A-34,  the  magnitude  of  D depends  on  the  distribution 
of  negative  ions.  Negative  ion  profiles  are  shown  in  Figure  2-9.  Note 
that  the  primary  reason  for  the  rapid  decrease  in  D for  altitudes  below 
about  60  km  is  the  shift  in  dominant  ions  from  0^  and  0 to  other  ions 
with  higher  electron  affinity.  The  recombination  rates  and  ou  are 
also  described  in  Appendix  A.  The  electron- ion  rate  depends  on  the 
distribution  of  positive  ions  (Equation  A-33) . Positive  ion  profiles 
arc  shown  in  Figure  2-10.  Note  that  the  primary  reason  for  the  increase 
in  coming  down  from  120  km  to  about  60  km  is  the  shift  in  dominant 
ions  from  those  with  slow  recombination  rates  to  those  with  higher  rates. 
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Figure  2-8.  Mean  daytime  effective  lumped  parameter  reaction  rate  coefficients 


MID-DAY 


CONCENTRATION  (cm'J) 

Figure  2-9.  Mean  daytime  natural  negative  ion  profiles. 


Ion  profiles. 


The  following  list  presents  the  ions  and  rates  in  order  of  in- 


creasing a^: 

Ion 

3 -1 

Recombination  rate  (cm  s ) 

(0*) 

k162 

= 2.  lx  10"7(T/300)  "°* 7 

(N0+) 

k157 

= 4. Ox  10"7(T/300)"1,0 

(0*) -group 

k164 

= 1.5 x 10~6(T/300) ”°* 7 

(O^fW)) -group 

k163 

='  : .5x  1G”6(T/300) “°* 7 

(N0+(N2)) -group 

k160 

= 1.5 x 10'6(T/300)'1,0 

(N0+(W)2) 

k158 

= 3.0 x 10"6(T/300)"0,2 

(N0+(W)3) 

k159 

= 5.0 x 10'6(T/300)'°‘2 

-6 

(»3o(W)2) 

k154 

= 5. lx  10 

-6 

(H30) (W)3) 

k155 

= 6.1 x 10  ° 

+ 

_ _ -6 

(H3o  (W)4) 

k156 

= 7.4  x 10 

The  lumped-parameter 

solutions  (Equations  A-36  through  A-39)  give  some 

insight  into  the  important  processes  in  the  absorbing  region.  From 
Figures  2-1  and  2-8,  note  that  D is  always  much  bigger  than  ^qa  . This 
means  that  in  Equation  A-57,  the  first  term  in  the  numerator  can  be 
neglected,  as  well  as  the  third  term  in  the  denominator.  Thus  the  ex- 
pression for  the  electron  concentration  reduces  to 

Co)  - 5Ta  <m*>t  (2-S) 

where 

(M+)t  =/q7 a . (2-6) 

Reflection  on  Figures  2-6,  2-7,  and  2-8  shows  that  the  crossover  be- 
tween D and  A determines  the  lower  boundary  of  the  absorbing  region 
and  that  through  most  of  the  absorbing  region  (e)  = yq/oij.  Thus  in 
the  sensitivity  analyses,  we  should  be  alert  to  processes  that  affect 
these  parameters. 


t 
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SECTION  2 


NUCLEAR  DISTURBED  IONOSPHERE 

Examples  of  enhanced  ion  production  rate  due  to  a nuclear  source 
are  shown  in  Figure  2-11.  The  nuclear  source  is  gamma-ray  energy  de- 
position from  radioactive  fission  debris.  The  sequence  of  cui'ves 
(those  labeled  A,  B,  C)  represent  the  relaxation  of  the  disturbance 
due  to  the  time  decay  of  the  radioactive  fission  products.  Each  curve 
represents  a factor  of  10  reduction  in  nuclear  energy  deposition  from 

the  previous  curve.  Since  the  decay  of  the  gamma-ray  energy  output 
-1  2 

follows  i (1 + t)  ‘ behavior  (t  is  time  after  burst),  the  sequence 
of  curves  represent  factors  of  about  7 increase  in  time  after  burst. 
These  curves  could  represent  relatively  fresh  fission  debris  from  a 
small  yield  burst  or  one  which  is  very  distant:  in  this  case  the  time 
intervals  between  curves  A,  B,  and  C would  be  relatively  short  (eg, 

A = 9 s,  B - 1 min,  C - 7 min).  The  curves  could  also  be  matched  by 
much  older  fission  debris  from  a larger  burst  or  one  nearby:  in  this 
case  the  time  interval  between  curves  could  be  much  longer  (eg,  A = 7 
min,  B = 50  min,  C = 6 hours).  The  curve  labeled  NATURAL  is  the  same 
as  that  in  Figure  2-1. 

The  corresponding  enhancement  of  electron  concentrations  and  HF 
wave  absorption  are  shown  in  Figures  2-12  and  2-13.  The  important 
thing  to  notice  in  Figures  2-11  through  2-13  is  that  the  nuclear  dis- 
turbance only  affects  the  ionosphere  below  about  80  km  altitude.  This 
result  is  representative  of  levels  of  disturbance  encountered  during 
the  recovery  phase  of  an  HF  system.  The  tic  marks  on  the  curves  in 
Figure  2-12  delineate  the  altitude  interval  where  90  percent  of  the 
total  vertical  integrated  absorption  occurs. 

From  Equations  2-5  and  2-6,  the  electron  and  total  positive  ion 
density  (and  therefore  also  the  total  negative  ion  density)  are  ap- 
proximately proportional  to  the  square  root  of  the  ion  production  rate 
q.  Variation  from  this  dependence  would  arise  only  from  a change  in 
the  effective  reaction  rates  A,  D,  a,  or  a.  as  a function  of  q.  Com- 
parison of  Figures  2-10  and  2-11  shows  that  the  square-root  dependence 
holds  to  a close  approximation.  Careful  examination  of  the  calculated 
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Enhanced  ion  production  rate  profiles  due  to  radio 
active  nuclear  fission  sources. 


CONCENTRATION  (cm  J) 

Figure  2-12.  Enhanced  electron  concentration  profiles 
radioactive  nuclear  fission  sources. 
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output  shows  that  in  the  nuclear  cases  D changes  at  most  about  1 per- 
cent from  that  shown  in  Figure  2-8,  due  to  slight  changes  in  the  rela- 
tive proportions  of  negative  ions  from  those  shown  in  Figure  2-9. 

(The  absolute  magnitude  of  each  ion  concentration  increases  in  the 

nuclear  cases,  of  course.)  The  value  of  the  electron-ion  recombination 

i 

rate  a changes  a little  more.  In  the  larger  (Case  A)  nuclear  case, 

> is  reduced  about  10  percent  from  that  shown  in  Figure  2-8  in  the 

altitude  regime  between  50  and  70  km.  This  is  due  to  slightly  higher 

proportions  of  ions  such  as  NO+(N~)  and  N0+(lV)  which  have  lower  rc- 

+ ^ + 

combination  rates  than  the  dominant  ions  11^0  (W) 7 and  H^O  (W)^.  The 
attachment  rate  A and  the  ion-ion  recombination  rate  cr  are  independent 
of  q and  remain  as  shown  in  Figure  2-8. 

The  simplified  interpretation  of  the  lumped  parameter  solution  for 
the  undisturbed  case  (see  Equations  2-5  and  2-6  above)  can  be  seen  to 
remain  valid  for  the  nuclear  disturbed  cases  also.  Even  with  the  larger 
values  for  q,  we  still  have  D » ^eja  throughout  the  absorbing  region. 
Thus,  the  crossover  between  D and  A determines  the  lower  boundary  of 
the  absorbing  region,  and  throughout  most  of  the  absorbing  region  we 
have  .(e)  = ^q/otj  . Thus,  in  the  sensitivity  analyses,  we  shall  address 
the  effect  of  the  minor  neutral  species  on  the  parameters  A,  D,  q,  and 


SECTION  3 

SENSITIVITY  ANALYSIS 


This  section  presents  an  investigation  into  the  sensitivity  of  the 
HF  wave  absorption  to  variations  in  the  atmospheric  minor  neutral  spe- 
cies, and  an  analysis  of  the  observed  sensitivity  in  terms  of  the  iono- 
spheric chemistry  model  parameters.  In  anticipation  of  this  investiga- 
tion, atmospheric  data  on  the  expected  variations  in  minor  neutral 
species  has  been  collected  and  collated  (Reference  3-1).  Unfortunately, 
most  of  the  available  data  is  for  the  stratosphere  (below  about  50  km 
altitude)  or  the  lower  thermosphere  (above  about  80  km  altitude).  In 
the  important  mesosphere  region  (50  to  80  km  altitude)  where  much  of 
the  HF  wave  absorption  takes  place — particularly  in  the  nuclear  disturbed 
cases — there  is  very  little  empirical  data.  Theoretical  models  for  the 
minor  neutral  species  concentrations  are  not  reliable  for  absolute  pre- 
dictions because  they  are  highly  sensitive  to  transport  processes  (molec- 
ulat  or  turbulent  diffusion,’  gross  wind  patterns)  that  are  difficult  to 
model.  Thus  we  have  chosen  somewhat  crude  bounds  for  the  minor  neutral 
species  concentration  variations  to  be  employed  in  the  sensitivity  analy- 
sis, based  on  a combination  of  extrapolated  empirical  data  and  reasoning 
from  the  theoretical  models. 

Since  ionospheric  absorption  of  HF  waves  is  significant  only  in  the 
daytime,  we  shall  concentrate  on  the  daytime  conditions.  Figure  3-1 
shows  mean  daytime  profiles  (Reference  3-2)  of  the  minor  atmospheric 
species  considered  in  the  analysis. 

Reference  3-1.  Jordano,  R.J.,  and  H.H.  Rutherford,  Suggested  Natural 
Variations  in  Atmospheric  Minor  Neutral  Species , GE78TMP-26  (ARBRL- 
CR-00375) , General  electric  Company— TEMPO,  July  1978. 

Reference  3-2.  Myers,  B.F.,  The  ROSCOE  Manual , Volume  14b -Midlatitude 
Density  Profiles  of  Selected  Atmospheric  Species , SAI-75-609-U-2B 
(DNA  3964l-'-14b) , Science  Applications,  Inc.,  June  1975. 
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Figure  3-1.  Mean  daytime  atmospheric  minor  neutral  species  profiles 
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OXYGEN  SPECIES  VARIATION 

The  minor  oxygen  species  of  concern  are  atomic  oxygen  (0) , o^one 
(0,) , and  excited  molecular  oxygen  in  the  a-state  ((^(a^Ag)).  The 
source  of  these  species  in  the  atmosphere  is  solar  photodissociation. 
The  principal  chain  of  reactions  are  as  follows.  Solar  photodissocia- 
tion of  molecular  oxygen  occurs  primarily  in  the  thermosphere  and  up- 
per mesosphere: 

C>2  + hv  -*•  0 + 0 (3-1) 

The  atomic  oxygen  then  recombines  with  C>2  to  form  0^  through  the  fol- 
lowing process: 

0 + 02+M-*-03  + M (3-2) 

where  M is  any  third  party.  Solar  photodissociation  of  0^  (this  ra- 
diation penetrates  down  into  the  upper  stratosphere)  then  produces 
the  excited  oxygen  molecule: 

Oj  + hv  -*■  02(a1Ag)  + 0 (3-3) 

Theoretical  treatments  of  diurnal  and  seasonal  variations  at  high 
latitude  may  be  found  in  References  3-3  and  3-4.  Below  about  80  km 
altitude,  the  (0)  and  (02(a1Ag))  disappear  at  night.  The  total  amount 
of  odd  oxygen,  (0)  + (0^)  , results  from  long-term  buildups  over  sev- 
eral diurnal  periods  and  depends  strongly  on  the  various  transport 
processes.  The  ratio  (0)/(C>3)  depends  strongly  on  the  balance  between 
Reactions  3-2  and  3-3,  with  the  altitude  dependence  on  the  three-body 
Reaction  3-2  accounting  in  large  measure  for  the  shapes  of  the  daytime 
concentration  profiles  in  the  mesosphere.  Since  the  Hartley  band 

Reference  3-3.  Shimazaki , T.,  and  A.R.  Laird,  "Seasonal  Effects  on 
Distributions  of  Minor  Neutral  Constituents  in  the  Mesosphere  and 
Lower  Thermosphere,"  Radio  Science  7_,  23(1972). 


Reference  3-4.  Koshelev,  V.V.,  "Diurnal  and  Seasonal  Variations  of 
Oxygen,  Hydrogen,  and  Nitrogen  Components  at  Heights  of  Mesosphere 
and  Lower  Thermosphere,"  J.  Atmo3.  Terr.  Phya.  38^  991(1976). 
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radiation  that  dissociates  the  0^  penetrates  below  60  km,  we  expect 


little  dependence  on  solar  zenith  angle  in  the  absorbing  region  except 
very  near  sunrise  or  sunset.  Thus  the  ratio  (0)/(03)  and  (O^a1  g))/ 
(Oj)  are  relatively  constant  throughout  the  day  and  there  should  be 
a positive  correlation  between  the  variations  in  these  species  (as 
seen  on  the  longer  time  scale  or  over  spatial  coordinates). 

Based  on  the  empirical  data,  we  have  chosen  a simple  variation  of 
plus  or  minus  a factor  of  3 (over  the  whole  profile)  to  be  representa- 
tive of  the  expected  latitude,  season,  or  solar  cycle  variations  in 
the  daytime  concentrations  of  the  minor  oxygen  species.  The  general 
trend  is  for  the  concentrations  to  be  reduced  at  higher  latitudes. 

There  is  some  confusion  as  to  the  seasonal  variations,  both  in  theory 
(see  References  3-3  and  3-4)  and  in  the  measurements  (see  Reference 
3-1). 

Table  3-1  presents  the  sensitivity  of  the  integrated  one-way  verti- 
cal HF  absorption  to  variations  in  the  oxygen  species.  Several  cases 


Table  3-1.  Integrated  one-way  vertical  10  MHz  absorption 
with  variation  of  oxygen  species. 


Species  Variation 

Natural  (dB) 

Nuclear  (dB) 

Base  Case 

1.41 

5.10 

3X  Increase  (03),  (02) 

1.43 

5.32 

3X  Increase  (0) 

1.60 

6.92 

3X  Increase  (0),  (0^) » (02) 

1.61 

7.04 

1/3 X Decrease  (03),  (02) 

1.40 

4.91 

1/3 X Decrease  (0) 

1.31 

3.86 

1/3 X Decrease  (0),  (03),  (0*) 

1.29 

3.62 

Note: 

(0*)  = (O^a^g)) 
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are  listed  which  are  variations  from  the  natural  ionosphere  case  shown 
in  Figure  2-6,  and  from  the  largest  nuclear  case  (Case  A)  of  Figure 
2-13.  In  each  case,  we  investigate  a variation  of  the  (0)  profile, 
the  (0j)  and  ((^(a^Ag))  profiles,  and  then  all  three  profiles  at  once. 
Note  that  the  sensitivity  due  to  the  (0)  variation  is  stronger  than 
that  due  to  (0-)  and  (0o(aaAg)),  and  that  the  sensitivity  is  larger 
in  the  nuclear  case  than  under  natural  conditions.  We  shall  analyze 
this  behavior  in  terms  of  the  model  parameters. 

Figure  3-2  shows  vertical  profiles  of  the  incremental  absorption 
for  the  natural  and  nuclear  cases  of  maximum  variation  in  Table  3-1. 

Note  that  most  of  the  sensitivity  occurs  at  altitudes  below  65-70  km, 
with  the  peak  variation  in  incremental  absorption  maximizing  at  about 
a factor  of  10  at  50  km  altitude.  This  altitude  regime  contributes  a 
small  portion  of  the  total  integrated  absorption  in  the  natural  case, 
but  contributes  a major  portion  of  the  integrated  absorption  in  the 
nuclear  case.  This  explains  the  larger  sensitivity  of  the  nuclear 
case  noted  in  Table  3-1. 

In  analyzing  the  results  further,  we  shall  employ  the  effective 
lumped  parameter  concepts  discussed  in  Section  2.  Recall  that  the  in- 
cremental absorption  is  proportional  to  the  electron  concentration  (e) , 
and  that 

<*>  - D^A  (3-4) 

through  most  of  our  regime  of  interest. 

In  the  altitude  regime  below  65-70  km,  the  oxygen  species  varia- 
tions act  through  the  effective  detachment  rate  D.  The  effects  of 
changes  in  D arc  washed  out  at  higher  altitudes  because  for  small  A 
the  ratio  D/(D+A)  goes  to  unity  in  all  cases.  Changes  in  D can  result 
from  changes  in  the' detachment  rates  for  particular  negative  ions,  or 
from  changes  in  the  relative  proportions  of  ions  with  different  detach- 
ment rates  (sec  liquation  A-34) . Figure  3-3  shows  the  contribution  of 
each  negative  ion  rate  to  the  total  effective  D for  the  base  case  condi- 
tion. (Figure  3-3  is  for  the  natural  ionosphere  case,  but  is  very 
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Figure  3-3.  Negative  ion  contributions  to  effective  detachment 
rate  D,  daytime  profiles. 
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representative  of  the  nuclear  cases  also  since  the  relative  distribu- 
tion of  negative  ions  is  a weak  function  of  q.)  Note  that  the  impor- 
tant ions  are  0~  and  0"  at  altitudes  where  the  HF  absorption  occurs. 
Other  ions  have  larger  concentrations,  but  have  negligible  detachment 
rates  (Equations  A-12  to  A-21). 

Figure  3-4  shows  the  contributions  of  specific  chemical  reactions 
to  the  detachment  rates  for  Ol  and  for  0”.  Note  the  strong  dependence 

l ^ 

on  (0)  and  (0-(a  Ag)).  Because  most  of  the  sensitivity  shown  in  Figure 

^ 1 
3-2  occurs  below  60  km,  we  would  expe  t the  (02(a  Ag))  variations  to 

have  more  impact  than  the  (0)  variations,  contrary  to  the  results  of 

Table  3-1.  The  answer  lies  in  the  different  ways  these  species  affect 

the  relative  distribution  of  negative  ions.  The  increase  of  (0_)  and 

1 - 
(02(a  Ag))  tend  to  decrease  the  relative  concentrations  of  0^  and  0 , 

compensating  for  the  increase  in  the  detachment  rates  of  those  species. 

The  increase  of  (0)  tends  to  increase  the  relative  concentrations  of 

0"  and  0~,  however,  enhancing  the  effectiveness  of  their  detachment 

1 

rates  in  contributing  to  D.  For  example,  increasing  (0^)  and  (02(a  g) 
by  a factor  of  3 increases  the  0~  and  0 detachment  rates  about  a fac- 
tor  of  2 to  3 between  50  and  60  km  altitude  (Figure  3-4) . This  also  • 
decreases  the  relative  concentrations  of  02  and  0 about  a factor  of 
2,  however  (base  case  distributions  are  shown  in  Figure  2-9),  leaving 
a net  increase  of  D of,  at  most,  50  percent.  Increasing  (0)  by  a fac- 
tor of  3 increases  the  0~  and  0~  detachment  rates  less  than  a factor 
of  2 between  50  and  60  km  altitude.  This  also  increases  the  relative 
concentrations  of  C>2  and  0 in  this  regime  by  a factor  of  1.5  to  3, 
giving  typical  increases  in  D more  than  a factor  of  3. 

It  is  difficult  to  identify  the  key  chemical  reactions  responsible 
for  these  shifts  in  relative  concentrations  of  02  and  0 , because  all 
the  ions  are  coupled  by  a network  of  transfer  reactions  (see  Figure 
A-2) . A survey  of  the  transfer  rates  involving  the  minor  neutral  spe- 
cies shows  the  following  candidates  for  the  production  of  C>2  and  0 : 
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DETACHMENT  RATE  (s_1) 

Figure  3-4.  Daytime  detachment  rate  profiles  for  (0")  and  (01). 


°3  * 0 - °2  * °2 

(k234  « 2.SX10"10) 

(3-5) 

COj  ♦ 0 0'2  ♦ C02 

(k203  * l.lxlO'10) 

(3-6) 

0“  ♦ 0 -»■  0"  ♦ 02 

(k22J  * 1.5  x 10'10) 

(3-7) 

Since  C02  is  a dominant  ion  in  this  regime,  the  second  reaction  (^203) 
is  most  likely  the  key  one.  The  third  reaction  (k223)  would  share  the 
increase  with  (0  ).  Two  candidate  re  ctions  contribute  to  the  loss 
of  0~  and  0“: 

°2  * °3  * °3  * °2 

(k226  - 4.0  «10-10) 

(3-8) 

°'  * °3  - °3  * 0 

(k218  * s.s.io-10) 

(3-9) 

The  first  reaction  (k^^)  contributes  30  to  40  percent  of  the  (0~) 
loss  rate  in  this  regime.  The  second  (k^g)  contributes  15  to  20  per- 
cent of  the  loss  rate  of  (0~) . 

We  make  one  final  note  with  regard  to  the  small  variations  in  in- 
cremental absorption  at  the  higher  altitudes,  above  90  km,  in  Figure 
3-2.  These  are  due  to  two  factors  related  to  the  atomic  oxygen  con- 
centration. One  is  the  direct  ionization  of  (0)  by  solar  UV  (Equation 
B-24) . This  represents  about  5 percent  of  the  total  UV  ionization 
source  at  above  100  km,  for  example  (see  Figure  2-1).  The  second  fac- 
tor is  the  influence  of  (0)  on  the  charge  transfer  from  (N*)  to  (0*) 

+.  z z 

or  (NO  ) (Equations  A-l  and  A-2,  specifically  terms  P and  Pfi)  . This 
has  negligible  effect  on  the  ratio  (NO  )/(0  ),  however,  because  the 
primary  source  of  (N0+)  is  from  charge  transfer  from  ot  in  cither  case. 

4.  + ^ 

(A  change  in  the  ratio  (NO  )/(02)  would  have  affected  a^) . 

ODD  NITROGEN  SPECIES  VARIATION 

The  minor  nitrogen  species  of  concern  are  atomic  oxygen  (N) , nitric 
oxide  (NO),  and  nitrogen  dioxide  (N02) . Figure  3-1  shows  mean  daytime 
profiles  for  these  species  concentrations  (Reference  3-2).  Potential 


— — -j- 
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sources  of  these  species  in  the  upper  atmosphere  are  ion-transfer  and 
ion-recombination  reactions,  and  solar  photodissociation  of  N2,  al- 
though the  quantitative  aspcccs  of  the  various  production  mechanisms 
are  not  well  understood  (Reference  3-5) . The  ion  sources  would  be  in 
the  thermosphere  (above  100  km  altitude)  and  would  include  the  follow- 
ing reactions: 


N*  ♦ e -*•  N(2D)  + N(4S) 


(3-10) 


N*  ♦ 0 -*■  NO*  ♦ N(4S  or  2D) 


0*  ♦ N2  -*•  NO*  ♦ N(4S  or  2D) 


NO*  ♦ e -►  0 ♦ N(4S  or  2D) 


(3-11) 


(3-12) 


(3-13) 


The  excited  atomic  nitrogen  then  reacts  rapidly  with  molecular  oxygen 
to  form  nitric  oxide: 


N(  D)  + 02  -*•  NO  + 0 


(3-14) 


The  principal  sinks  for  odd  nitrogen  at  these  altitudes  are  transport 

8 2 

processes  (downward  flux  of  (NO)  into  the  mesosphere  exceeds  10  cm 
s"1)  and  the  following  chemical  reaction 


(N  S)  + NO  -*•  N2  + 0 


(3-15) 


From  the  last  two  reactions,  we  see  that  the  relative  abundance  of 
ionospheric  (N)  versus  (NO)  depends  critically  on  the  assumed  branch- 
ing ratios  (to  N(4S)  or  N(2D))  in  the  ion  reactions.  Most  theoretical 
2 

models  favor  the  (N(  D)),  and  thus  give  very  low  ionospheric  concen- 
trations of  (N)  compared  to  (NO)  below  120  km  altitude.  F.inpirical 


Reference  3-5.  Strobcl,  O.F.,  "Minor  Neutral  Constituents  in  the 
Mesosphere  and  Lower  Thermosphere,"  Radio  Science  7j  1 (1972). 
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measurements  of  (N)  have  not  substantiated  this  condition.  In  fact, 
they  give  (N)  concentrations  typically  3 orders  of  magnitude  higher 
than  the  calculations  (Reference  3-1).  The  measurements  are  difficult 
however,  because  the  in  aitu  techniques  arc  subject  to  contamination. 

In  the  mesosphere,  the  important  reactions  (assuming  (N)  is  negli- 
gible) are: 


NO  ♦ 0, 


no2  ♦ o2 


N02  + 0 -*•  NO  ♦ 02 


(3-16) 


(3-17) 


N02  + hv  NO  + 0 


(3-18) 


NO  + 0 ■*>  N02  + hv 


(3-19) 


NO  + 0 + M -*•  N02  + M 


(3-20) 


Note  that  all  the  reactions  just  convert  between  (NO)  and  (N02)  so 
that  (on  a time  scale  shorter  than  the  diurnal)  the  total  (NO)  * (N02) 
is  constant.  The  ratio  (N0)/(N02)  is  controlled  by  solar  photodis- 
sociation of  (N02)  and  by  the  ratio  (0)/(03).  Below  70  km,  the  NO  is 
converted  to  N02  at  night.  Theoretical  calculations  (Reference  3-3) 
show  that  above  about  50  km  altitude,  the  (NO)  is  fairly  constant 
throughout  the  day,  with  the  diurnal  transition  occurring  abruptly  at 
sunrise  and  sunset.  The  (N02)  shows  a stronger  dependence  on  solar 
zenith  angle  during  the  day  due  to  the  photodissociation  process. 

Thus,  time  of  day  must  be  considered  when  comparing  latitudinal,  sea- 
sonal,  or  solar  cycle  variations  in  (N02) . For  a given  solar  zenith, 
however,  we  expect  positive  correlation  between  the  variations  in  (NO) 
and  (N02) . 

Empirical  data  (Reference  3-1)  indicate  that  variation  in  (NO)  and 
(N02)  are  probably  larger  than  those  in  (0)  and  (0j) . There  is  a de- 
finite increase  in  (NO)  and  (N02)  at  high  latitudes,  which  appears  to 
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be  related  to  auroral  activity  (Reference  3-6) . Seasonal  data  tend 
to  favor  maximum  concentrations  in  summer.  To  simplify  comparison  to 
the  oxygen  snccics  sensitivity,  we  have  varied  the  (NO)  and  (NO^)  by 
the  same  factor  of  3 in  this  portion  of  the  analysis.  We  shall  dis- 
cuss the  potential  effect  of  the  very  large  uncertainty  in  (N)  toward 
the  end  of  this  section. 

Table  3-2  presents  the  sensitivity  of  the  integrated  one-way  ver- 
tical I1F  absorption  to  variations  in  the  (NO)  and  (NO2) . The  cases 


Table  3-2.  Integrated  one-way  vertical  10  MHz  absorption 
with  variation  of  nitrogen  species. 


Species  Variation 

Natural 

(dB) 

Nuclear 

(dB) 

Base  Case 

1.41 

5.10 

3 X increase  (NO),  (N02) 

2.05 

5.52 

1/3 X decrease  (NO),  (N02) 

1.10 

4.91 

listed  are  variations  of  the  natural  ionosphere  case  shown  in  Figure 
2-6,  and  from  the  largest  nuclear  case  (Case  A)  of  Figure  2-13.  In 
contrast  to  the  observed  sensitivity  to  oxygen  species  (Table  3-1) , 
the  effect  of  the  nitrogen  species  variation  is  stronger  for  the  natu- 
ral ionosphere  than  for  the  nuclear-disturbed  ionosphere. 

Figure  3-5  shows  vertical  profiles  of  the  incremental  absorption 
for  each  of  the  cases  of  Table  3-2.  Note  that  the  largest  sensitivity 
is  in  the  altitude  band  between  65  and  95  km  and  there  is  a crossover 
in  the  correlation  above  100  km  altitude.  In  terms  of  the  lumped  param- 
eters (liquation  3-4),  the  impact  of  the  nitrogen  species  is  principally 
through  the  ion  production  rate  q.  There  is  some  smaller  effect  through 
the  electron  recombination  rate  ctj,  and  an  almost  negligible  effect 
through  the  detachment  rate  D. 

Reference’  3-6. ~ Gerard,  .J.C.,  and  C.A.  Barth,  "Hi gh-I.at i tilde  Nitric 

Oxide  in  the  l.ower  Thermosphere,"  J.  Ceophj3.  /fes.  #2,  674  (1977). 
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The  effect  on  the  production  rate  q is  entirely  through  the  solar 
Lyman-a  ionization  (Equation  B-4)  which  is  directly  proportional  to 
(NO).  Figure  2-1  shows  that  Lyman-a  dominates  the  natural  ionization 
sources  between  65  and  95  km  altitude.  Figure  2-10  shows  that  the 
nuclear  ionization  (Case  A)  would  tend  to  mask  sensitivity  to  Lyman- 
below  about  75  km,  however,  where  most  of  the  HF  wave  absorption  takes 
place  in  the  nuclear  case.  This  explains  the  reduced  sensitivity  to 
the  nitrogen  species  variations  in  the  nuclear  case. 

The  effect  of  the  nitrogen  species  variation  on  the  electron  re- 
combination rate  a,  is  through  the  positive  ion  transfer  coefficient 

Q + + 

Tp^  (Table  A-l),  which  enhances  the  production  of  NO  from  0^  through 

the  following  reactions: 


+.  N -*>  N0+  +0 

<k13S  * 12 * 10'10> 

(3-21) 

+ NO  ->  N0+  + 02 

(kU6  » “.S*!0’10) 

(3-22) 

+ \'2  ->  N0+  + NO 

(k139  = 1.0  xio'16) 

(3-23) 

With  the  base  case  distribution  of  species  (Figure  3-1),  the  second 
reaction,  that  involving  (NO),  dominates  above  about  85  km,  and  the 
third  reaction,  with  (N2)  , is  most  important  below  about  85  km.  Above 
90  km,  where  the  (NO)  variation  has  the  strongest  effect,  the  dominant 
ions  are  0*  and  N0+  (see  Figure  2-10).  With  increased  (NO),  the  en- 
hanced transfer  from  0*  to  N0+  increases  a^  about  10-20  percent  (since 
the  N0+  recombination  rate  is  larger  than  that  for  0*) . The  electron 
density  is  thus  decreased  about  5-10  percent. 

The  variation  of  nitrogen  species  has  some  effect  on  the  detach- 
ment rate  D.  Increasing  (NO)  increases  the  0~  and  0~  rates  slightly 
(Equations  A-12  and  A-13) , but  the  major  effect  is  the  shifting  of  the 
relative  negative  ion  concentrations  to  favor  N02  , NO^  , and  00N0  . 
The  increased  (NO)  and  (N02)  has  a net  effect  of  a 10-20  percent  de- 
crease in  D in  the  altitude  band  between  50  and  65  km  where  the  value 
of  D is  important. 
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We  now  address  the  sensitivity  to  the  atomic  nitrogen  (N) . We 
have  not  displayed  the  results,  but  the  effect  of  (N)  on  the  HF  ab- 
sorption is  entirely  through  the  positive  ion  transfer  of  Equation 
3-21.  The  rate  coefficient  is  about  one-fourth  that  for  the  transfer 
involving  (NO)  (Equation  3-22)  and  one  can  see  from  Figure  3-1  that 
increasing  (N)  about  two  orders  of  magnitude  would  make  it  competitive 
with  (NO) . The  net  effect  of  the  very  large  (N)  would  then  be  to  shift 
the  relative  ion  concentrations  from  0*  to  N0+  above  about  95  km  alti- 
tude, increasing  and  decreasing  the  electron  concentration  and  H- 
wave  absorption  in  this  region.  The  maximum  possible  reduction  in 
absorption  due  to  this  effect  would  be  about  10-15  percent  reduction 
of  that  portion  above  95  km.  This  is  an  insignificant  effect  even  in 
the  natural  ionosphere  case. 

WATER  VAPOR  .AND  CARBON  DIOXIDE  VARIATIONS 

There  are  insufficient  data  upon  which  to  base  estimates  of  expec- 
ted variations  in  thes  species  (Reference  3-1) . In  the  mesosphere  and 
lower  thermosphere,  where  most  HF  wave  absorption  occurs,  there  are 
almost  no  empirical  measurements  and  estimates  of  the  concentrations 
are  based  on  model  calculations.  Consequently,  we  have  chosen  to  use 
a simple  factor  of  3 variation  of  the  profiles  of  Figure  3-1  to  facili- 
tate comparison  of  sensitivities  to  the  other  species. 

Table  3-3  presents  the  sensitivity  of  the  integrated  one-way  verti- 
cal mp  absorption  to  variations  in  the  (H20)  and  (CO,,)  profiles.  The 
listed  are  variations  of  the  natural  ionosphere  case  shown  in 
Fig  Te  2-6,  and  from  the  largest  nuclear  case  (Case  A)  of  Figure  2-13. 
Note  that  increases  in  either  the  (H20)  or  (CO.,)  cause  a decrease  in 
the  HF  absorption.  Both  the  (H.,0)  and  (C02)  variations  cause  about 
5 percent  variation  in  absorption  in  the  natural  ionosphere  case.  In 
the  nuclear  case,  the  effect  of  the  (C02)  variation  is  reduced  slightly, 
and  the  effect  of  the  (H20)  variation  is  enhanced  slightly. 

Figure  3-6  shows  vertical  profiles  of  the  incremental  absorption 
for  each  of  these  cases.  Note  that  the  sensitivity  to  (C02)  is  almost 
entirely  in  the  short  altitude  band  between  75  and  90  km,  while  the 
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absorption  profiles  with  variation  of  (H,0)  and  (C(L). 
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Table  3-3.  Integrated  one-way  vertical  10  MHz  absorption 
with  variation  of  (I^O)  and  ( CO2 ) - 


Species  Variation 

Natural 

(dB) 

Nuclear 

(dB) 

Base  Case 

1.41 

5.10 

3X  increase  (H^O) 

1.36 

4.79 

1/3 X decrease  (H^O) 

1.49 

5.82 

3X  increase  (CO^) 

1.36 

4.88 

1/3 X decrease  ( CO2 ) 

1.48 

5.25 

sensitivity  to  (ll20)  variation  is  at  altitudes  between  55  and  80  km. 

The  absorption  in  the  nuclear  case  is  thus  more  sensitive  to  (I^O) 
because  the  nuclear  enhancement  is  entirely  within  the  lower  region. 

The  impact  of  the  (C02)  variation  is  through  (see  liquation  3-4). 
The  important  reaction  is 

NO+  + C02  + M -*  N0+(C02)  + M (3-24) 

Larger  (C0o)  enhances  the  transfer  from  (NO+)  to  (N0+(C02))  which  has 
a larger  electron  recombination  rate  (by  a factor  of  4,  compare  k^^ 
and  in  Appendix  A).  The  (NO+(CC>2))  is  lumped  with  the  (N0+(N2))- 

group  in  Figure  2-10.  They  are  lumped  together  for  convenience  in  the 
model  because  they  have  similar  recombination  rates  (see  liquation  A-30) . 
From  Figure  2-10  then,  we  see  that  the  important  interplay  between  the 
relative  concentration  of  N0+(CC>2)  and  other  ions  is  in  the  altitude 
regime  75-90  km. 

The  (C02)  also  has  influence  on  the  shift  in  negative  ion  concen- 
tration to  favor  (CCLj)  and  (CO^)  instead  of  (0  ) and  (C>2) . See  T^^ 

and  T „„  in  Tabic  A-2.  This  decreases  the  effective  detachment  D, 
n.3a 

but  the  magnitude  of  the  change  is  negligible  in  the  model. 
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The  impact  of  the  (H^O)  variation  is  also  through  a^.  Dozens  of 
positive  ion  transfer  coefficients  involving  (l^O)  favor  the  transi- 
tion to  higher  order  hydrates  of  the  ions  N0+  and  11^0*.  (See  Table 
A-l) . Higher  order  hydrated  ions  tend  to  have  larger  electron  recom- 
bination rates,  as  shown  in  Section  2 (see  list  prececding  Equation 
2-5) . The  largest  rates  are  more  than  an  order  of  magnitude  larger 
than  those  for  the  simple  ions  0*  and  N0+.  From  Figure  2-10,  we  see 
that  the  various  hydrated  ions  begin  to  become  important  below  about 
SO-km  altitude,  consistent  with  the  o..sct  of  the  (II20)  sensitivity  in 
Figure  3-6. 


APPENDIX  A 

THE  CHEMISTRY  MODEL 


The  chemistry  model  employed  for  calculations  in  this  report  is 
the  steady-state  ion  chemistry  model  taken  from  the  ROSCOE/WEPH  com- 
puter programs  (References  A-l  and  A-2) . We  summarize  the  properties 
of  the  model  here  for  the  reader's  convenience.  Note  in  particular 
that  neutral  species  are  inputs  for  this  part  of  the  model,  and  are 
assumed  to  be  unaffected  by  the  ion  solutions. 

The  positive  and  negative  ions  treated  explicitly  in  the  model  are 
illustrated  in  Figures  A-l  and  A-2,  along  with  a schematic  representa- 
tion of  the  positive  ion  and  negative  ion  transfer  reactions.  The 
numbers  next  to  the  arrowheads  indicate  the  ion  transfer  coefficients 
which  are  defined  in  Tables  A-l  and  A-2.  Additional  reactions  not 
shown  in  Figures  A-l  and  A-2  include  positive  ion  production  due  to 
natural  or  nuclear  energy  deposition,  the  production  of  negative  ions 
by  electron  attachment  to  neutral  species,  the  loss  of  negative  ions 
due  to  electron  detachment,  and  recombination  of  positive  ions  and 
electrons.  All  reactions  and  rate  coefficients  are  given  in  Table  A-3. 

The  ion  production  models  due  to  natural  or  nuclear  energy  deposi- 
tion provide  "primary"  ion  production  rates  q(N+),  q(0+),  qCN^),  qCO^), 
and  q(NO+) . The  first  three  species  (N+,  0+,  and  N*)  are  assumed  to 

Reference  A-l.  Knapp , W.S.,  et  al,  WEPH  VI:  A Fortran  Code  for  the 
Calculation  of  Ionization  and  Electromagnetic  Propagation  Effects 
Due  to  Nuclear  Detonations  (U)j  Volume  Z — Computational  Models  (U), 
GE75TMP-53  (DNA  3766T-3) , General  Electric  Company— TEMPO,  October 
1976  (Confidential/Formerly  Restricted  Data) . 

Reference  A-2.  Knapp,  W.S.,  et  al.  The  ROSCOE  Manual:  Volume  II — At- 
mospheric Chemistry  Models , GE74TMP-59  (DNA  3964F-11),  General 
Electric  Company — TEMPO  (to  be  published). 
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NOTE:  NUMBERS  REFER  TO  THE  POSITIVE  ION  SPECIES 


ANO  POSITIVE  ION  TRANSFER  COEFFICIENTS 

Figure  A-l.  Positive  ion  transfer  coefficients. 
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Figure  A-2.  Negative  ion  transfer  coefficients. 
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Figure  A-2.  Continued  (notes). 
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41 

13 

42 

43 

44 

14 

45 

* 

1 

°2 

6 

co- 

11 

CO’^O) 

2 

0* 

7 

NO^ 

12 

no2(h2o) 

3 

00N0' 

8 

NO3 

13 

no‘(h2o) 

4 

°; 

9 

CO4 

14 

co'(h2o) 

5 

°3 

10 

0^0) 

Table  A-l.  Positive  ion  transfer  coefficients. 


Tpi  -w«  + w°> 

Tp2  = ki47(02(lA)) 

Tp3  = k141<M> 

Tp4  = ki42^°2^M^ 

Tp5  = k150^V 
Tp6  = k148(02) 

Tp7  = ki34(H2°^M) 

Tp8  = k149^H20^ 

Tp9  = k109^ 

TplO  = k126^H20^ 

Tpll  = k106^H2°^M^ 

Tpl2  = ki08^H2°^M^ 

Tpl3  = klll^ 

Tp 14  = k 135  ^ ^ ^ + k136^N0^  + k139^N2  ^ + k138^N02^ 
Tpl5  = n151(N0) 

Tpl6  = WN0> 

Tpl7  = k118(M^ 

Tpl8  = k129(N0) 

Tp 19  = k115^ 

Tp20  = kH7(N2)(M^ 

Tp21  = k 125 

Tp22  = ki22^H2°^M^ 
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Table  A-l.  (Continued). 


Tp23 

= 

k124(H20)(M) 

Tp24 

* 

kno(H20)(H) 

Tp25 

k140(N2)(M) 

Tp26 

= 

k144(H20) 

Tp27 

= 

k105^M) 

Tp28 

= 

W> 

Tp29 

k103^ 

Tp30 

k12Q(HO)  + k121(H02) 

Tp31 

= 

k145(H2°) 

Tp32 

= 

k107^) 

Tp33 

= 

k100(H20)(M) 

Tp34 

k104^H20^ 

Tp35 

=2 

k105A^H20) 

Tp36 

k102(N2)(M) 

Tp37 

= 

ki37(N02) 

Tp38 

= 

kl52(N02) 

Tp39 

= 

kn4(C02)(M) 

Tp40 

= 

k 127 ( C02 ) 

Tp41 

= 

k116(H20)(M) 

Tp42 

= 

kl28(H2°) 

Tp43 

= 

k123(M) 

Tp44 

= 

kng(H20) 
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Table  A-2.  Negative  ion  transfer  coefficients 


k225^  + 

°sl  ,4 

^234(0) 

1^203(0) 

°sl,9 

Dsl,10 

k221(N)  + 

k223(0) 

°s2,6 

k238(N0) 

k206(N0) 

k228(N0) 

k224(02)(M) 

k236(°2) 

k226(03) 

k2 1 7 ( 02 ) ( M ) + k2j8(0|j) 

k239(0)  + 

k240(03) 

k208^°3^ 

k229(o3) 

k215(C02)(M) 

photodissociation  rate  producing  ith  negative 
ion  from  jth  negative  ion. 


Tn23 

= k222^®2^ 

Tn24 

= k216^N02^ 

Tn25 

= k241^^ 

Tn26 

= k231(N0)  + k232(N02*) 

Tn27 

k2oi(NO) 

00 

CM 

c 

»- 

* Ds7 ,12 

Tn29 

k233^NtV 

Tn30 

= k202^2^ 

Tn31 

= k211(N02)  + k212(03) 

^32 

= Ds8,13 

Tn33 

= k219(C02)(M) 

Tn34 

= k235(C02) 

Tn35 

= k227(C02) 

Tn36 

= Ds9,14 

Tn37 

= k22Q(H20)(M) 

Tn38 

= k237(H20) 

Tn39 

= k2Q0(H20)(M) 

Tn40 

= k21Q(N0) 

Tn41 

- k204(H20)(M) 

Tn42 

- k214(H20)(M) 

Tn43 

= k214^°3^ 

Tn44 

“ k2Q9(NO) 

Tn45 

= k205(H20)(M) 

' 


L i 

t 


I 


I 


f 

L 


I 


produce  0*  and  NO*  instantaneously  and  are  not  further  considered  ex- 
plicitly in  the  model.  The  effective  0*  and  NO*  production  rates  are 
computed  from 


qe(02)  * q(0*)  ♦ P^N*)  ♦ P3q(0+)  ♦ Psq(N*) 
qe(N0+)  = q(N0*)  + P2q(N+)  ♦ P^O*)  ♦ P6q(N*) 

where 

p - km 

1 k112  * k113 


(A-l) 

(A-2) 

(A- 3) 


1 - P1 

(A-4) 

k133(02) 

k131^V  * k132^M^N2^  * k133^V 

(A-r5) 

1 " P3 

(A-6) 

k130(°2) 

k130<V  * k180^ 

(A-7) 

i-^  . 

(A-8) 

Reactions  and  rate  coefficients  are  defined  in  Table  A-3. 


Electron  attachment  to  02  , 0^  , and  N02  produce  02  , 0",  and  NO, 


attachment 

rates 

(s'  ) are  computed  from 

(O')  : 

A1  - 

k166<N2><°2>  * k168(02)(°2> 

+ k179(02)(H20) 

(A-9) 

(O')  : 

> 

ro 

ii 

k170<°3>  * k178<0) 

(A-10) 

(N02)  : 

it 

to 

< 

kl65(N°2)  • 

(A-l 1) 

I 


1 
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Detachment  processes  include  collisional  and  associative  detachment 


from  O'  and  O' 

and  photodctachmcnt  from  Ol  , O',  0* 
- 2 2 3 

The  detachment 

rates  (s  ) are  computed  from 

(Oj) 

: D,  - k167(N2)  . k1M(02)  . k17500 

* k177<V'4»  * DU 

(O') 

: D2  = k1?1(N)  + k1?2(NO)  + k1?3(0) 

(OjJ 

: °5  = DtS 

(NOp 

: °7  - Dt7 

(NO') 

: °8  = Dt8. 

where  the  solar  photodetachment  terms  Dtl  are  taken 

(O’) 

: 0tl  - 0.3 

(O') 

: Dt2  = 1.4 

(03) 

: Dts  = 0.08 

(NOj) 

: Dt?  = 0.08 

(NO') 

: Dt8  = 0.03 

176 


(A-12 


(A-13 
(A- 14 


(A-1S 
(A- 16 


-1, 


(A-17) 
(A- 18] 
(A-19J 
(A-2o; 
(A-21; 


Very  little  data  is  available  on  ion-ion  recombination  coefficients, 
and  the  rate  is  assumed  independent  of  the  particular  positive  or  nega 
tivc  ion.  The  ion- ion  recombination  rate  coefficient  is  computed  from 


°i  * k242  + k243(M)  C",3s‘1 


(A-22] 
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The  electron-ion  recombination  coefficients  are  taken  as  follows 

, 3-1, 

(cm  s ) 


‘dl 

= k162 

(o2) 

ld2 

vO 

^■4 

M 

II 

(04 

ld3 

= k163 

(o20 

and  1 

!d4 

= k154 

0*3°' 

ld5 

= k155 

(h3° 

(A-23) 

(A-24) 


i30+(n2] 


(h3o  (W)4) 


(N0+) 


(A-25) 
(A- 26) 

(A-27) 

(A-28) 

(A-29) 


(NO+(N2),  N0+(C02),  and  NO+(W))  A-30) 


(NO  (W)2) 


(NO  (W)3) 


(A-31) 

(A-32) 


The  steady-state  solution  is  obtained  by  writing  the  differential 
equation  for  each  ion  species  and  setting  the  time  derivative  for  each 
species  concentration  to  zero.  This  results  in  a set  of  coupled,  non- 
linear equations  and  iterative  techniques  must  be  used  for  solution. 

The  procedure  employed  makes  use  of  the  fact  that  if  the  total  electron 
and  ion  concentrations  are  known,  the  distribution  of  ions  (the  ratio 
of  one  ion  concentration  to  another)  can  be  found  algebraically.  If 
the  ion  distributions  are  known,  effective  lumped  parameter  reaction 
rate  coefficients  can  be  defined  to  solve  for  the  total  ion  concentra- 
tions as  follows: 


llllhfrira  i ni  ii 
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“d  = 

(A-33) 

D = 

E D.(M')./(M')T 

(A-34) 

A = 

E A, 

(A- 35) 

1 

where  a, . are 
di 

the  various  electron- ion  recombination  rates 

(Equations 

A-23  through  A-32) , are  the  various  detachment  rates  (Equations  A-12 


through  A- 16),  A.  arc  the  various  attachment  rates  (Equations  A-9 
through  A-ll),  (M  ) ^ is  the  concentration  of  a particular,  positive 
ion,  (M  )T  is  the  total  positive  ion  concentration,  and  (M  )^  and  (M")t 


are  particular  and  total  negative  ion  concentrations,  respectively. 


The  effective  ion-ion  recombination  rate  or  has  been  assumed  independent 


of  the  particular  ions  (Equation  A-22) . With  the  effective  rates,  the 
lumped  parameter  solutions  are 


(M+).r  = v/qT^T 


(A-36) 


(e) 


q + D(M  ), 


A + D + ad(M  )T 


(A-37) 


(M')t  = (M  )T  - (e) 


(A-38) 


where  a is  obtained  from  the  cubic  equation 


Aa.  + 
1 


Dot, 


a. a. 
1 d 


v/q/a 


A + D + 


or  y/qTa 


(A-39) 


and  q is  the  total  ion  production  rate.  Thus,  in  the  iterative  pro- 
cedure, the  ion  concentration  ratios  are  calculated  from  the  total 
electron  and  ion  concentrations,  effective  lumped  parameter  reaction 
rates  arc  defined  from  the  ion  concentration  ratios,  and  the  total 
electron  and  ion  concentrations  are  computed  from  the  effective  rates. 
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The  iteration  converges  rapidly  because  the  ion  ratios  are  relatively 
insensitive  to  the  total  electron  and  ion  concentrations.  The  expres- 
sions for  the  ion  ratios  are  quite  complex  and  will  not  be  presented 
here. 

The  following  reactions  and  reaction  rate  coefficients  are  employed 
in  the  chemistry  model.  The  coefficients  A,  B,  and  C are  given  for  the 
reaction  rate  formula 


where  T is  the  air  temperature  (k) . 


Table  A-3.  Reactions  and  rate  coefficients. 


Reactions 

Rate  Coefficients 

ABC 

Neutral  Species 

Reactions 

1 

0^  + hv 

o2  + 0 

2 

0^  + hv 

->■ 

O^A)  + 0 

(Effective  values  com- 
puted within  model) 

3 

C>2  + hv 

0 + 0 

4 

N02  + hv 

NO  + 0 

* 5 

n(2d)  + o2 

-* 

NO  + 0 

7 . 5 ( - 12 ) +0.5 

6 

o2(1a) 

02  + hv 

2.6(-  4) 
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Table  A-3.  (Continued). 


Reactions 

Rate 

A 

Coefficients 

B 

C 

Neutral  Species  Reactions 

* 7 

o2(1a)  + 

N - 

02  + N 

2.0(-14) 

600 

* 8 

02(1A)  + 

N -*■ 
N2 

°2  + N2 

0 

* 9 

o2(1A)  + 

o -*■ 

o2  + 0 

0 

10 

O^A)  + 

°2 

°2  + °2 

2 . 2 ( - 18 ) 

+0.8 

11 

o^a)  + 

°3 

02  + 02  + 0 

4.5(-ll) 

2800 

*12 

N + NO 

-*> 

N2  + 0 

4.  l(-ll) 

410 

*13 

N + NO 

-► 

N2  + 0(lA) 

4.1(-11) 

410 

*14 

N + N02 

NO  + NO 

0 

*15 

N + N02 

-+■ 

n2o  + 0 

2 . 0( - 11 ) 

800 

16 

N + 0 

NO 

1 - 5( —17) 

17 

N + 0 + M 

NO  + M 

1 . 0( -32 ) 

-0.5 

18 

N + 02 

-► 

NO  + 0 

3. 3( -12) 

+1.0 

3150 

*19 

N +03 

NO  + 02 

3. l(-ll) 

1200 

20 

NO  + 0 

-► 

N0?  + hv 

6 . 6 ( - 17 ) 

-1.9 

*21 

NO  + 0 + 

M - 

N02  + M 

6 . 4 ( -32 ) 

-2.0 

22 

NO  + 03 

no2  + o2 

9.0(-13) 

1200 

23 

no2  + 0 

NO  + 02 

9 . 1 ( - 12 ) 

24 

o + o2  + 

M - 

<VM 

1 • 1 (-34) 

-510 

*25 

0 + 03 

-¥ 

O^A)  ♦ 02 

0 

*26 

0 + °3 

-¥ 

02*°2 

1.9(-li) 

2300 

*27 

o + ho2 

-► 

H + 02 

4.2(-ll) 
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Table  A-3.  (Continued). 


Reactions 


Rate  Coefficients 


Neutral  Species  Reactions 

*28 

o + ho2 

-*•  OH  + 02 

8 . 0 ( - 1 1 ) 

500 

*29 

o + h2 

- H + OH 

9 . 0( - 12 ) 

+1.0 

4500 

*30 

OH  + H2 

- H + H20 

3.6(-ll) 

2590 

*31 

OH  + N 

- H + NO 

5.3(-ll) 

*32 

H + 03 

- OH  + 02 

2 . 6 ( - 1 1 ) 

*33 

H + 03 

h-  ho2  + 0 

2 - 6 ( - 1 1 ) 

*34 

OH  + 03 

->  H02  + 02 

1 - 6 ( - 12 ) 

1000 

*35 

H + 02  + M 

- H02  + M 

5 . 5 ( -32 ) 

-1.0 

*36 

H + N02 

- OH  + NO 

5 . 8 ( - 10 ) 

740 

*37 

H02  + NO 

-*■  OH  + N02 

2 . 0( - 13 ) 

*38 

ho2  * 03 

-*■  OH  + 02  + 02 

1 .0( -13) 

1250 

*39 

NO  + OH  + M 

- HN02  + M 

6. 1 ( -31 ) 

-2.0 

*40 

N02  + OH  + M 

HN03  + M 

1.8C-30) 

-2.0 

*41 

no2  + ho2 

->  HN02  + 02 

3 . 0 ( - 14 ) 

*42 

OH  + H02 

-»■  H20  + 02 

2.0(-ll) 

*43 

o(!a)  + h2o 

- OH  + OH 

2 . 4 ( - 10) 

*44 

0(1A)  + N2 

->  0 + N2 

1 • 3 ( - 1 1 ) 

-260 

*45 

0( XA)  + 02 

0 + 02 

2 . 3 ( - 1 1 ) 

-170 

*46 

0(JA)  + H2 

-*  H + OH 

1 . 3 ( - 10) 

*47 

0( lA)  + H202 

-*•  OH  + H02 

2 . 5 ( - 10) 

*48 

H + H02 

> h2  + o2 

4 - 2 ( - 1 1 ) 

350 

*49 

H + H02 

- H20  + 0 

8. 3( -1 l) 

500 
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Table  A-3.  (Continued). 


Reactions 


Positive  Ion  and  Electron  Reactions 


*100  H30  +W  + M 
*101  H+(W)2  + M 

*102  H30++N2  + M 

103  H30+(N2)+M 


H+(W)2  + M 
h3o+  + W + M 
H30+(N2)  + M 
h3o+  + n2  + M 


Rate  Coefficients 


3 . 4 ( -27 ) -2.0 

1.6(+  1)  -2.0 
0 
0 


104 

h3o+(oh)  + W 

- H+(W)2  + OH 

1 -4(—  9) 

*105 

h3o+(oh)  + m 

- H30+  + OH  + M 

3.0(-  3) 

-2.0 

105A  H30+(N2)  +H20 

- h+(w)2  + N2 

0 

106 

h+(w)2  + w + m 

-*•  H+(W)3  + M 

2 . 3( -27 ) 

-2.0 

*107 

H+(W)3  + M 

- H+(W)2  + W + M 

1 - 6( - 1) 

-2.0 

108 

h+(w)3  + m + w 

- H+(W)4  + M 

2 . 4( -27 ) 

-2.0 

109 

H+(W)4  + M 

- H+(W)3  + W+M 

I.0(-  I) 

-2.0 

*110 

h+(w)4  + w + m 

- H+(W)5  + M 

9.0C-28) 

-2.0 

*111 

H+(W)5  + M 

- H+(W)4  + W + M 

1.3(+  0) 

-2.0 

112 

N+  + 02 

- 02  + N(20) 

2 . 8( - 10) 

113 

N+  + 02 

- N0+  + 0 

2.8(-10)#. 

*114 

N0+  + C02  + M 

■*  N0+(C02)  + M 

3 . 0( -29) 

-2.0 

*115 

N0+(C02)  + M 

- N0+  + C02  * M 

3.0(-  5)  '* 

-2.0 

116 

N0+  + W + M 

- N0+(W)  + M 

1 . 5 ( -28) 

-2.0 

*117 

N0+  + N2  + M 

N0+(N2)  + M 

2.0(-31) 

-4.4 

*118 

no+(n2)  + M 

- N0+  + N2  + M 

i.l(-  8) 

-4.4 

119 

no+(co2)  + W 

- N0+(W)  + C02 

o\ 

1 

o 

*120 

N0+(W)  ♦ OH 

- h3o+  + no2 

6.0(-ll) 

18200 


12000 


11300 
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Table  A-3.  (Continued). 


Rate  Coefficients 


*121  N0+(W)  + H02 

122  N0+(W)+WH4 

*123  N0+(W)2  + M 

124  N0+(W)2  + W + M 

125  N0+(W)3  <■  M 

126  N0+(W)3  + W 

127  N0+(N2)  + C02 

128  N0,'(N2)  + H20 

129  NO*  + NO 

130  N2  + 02 

131  0*  + N2 

*132  0*  + N2  + M 

133  0*  + 02 

134  0*  + W + M 

135  02  *■  N 

136  02  t-  NO 

137  02  > N02 
*138  02  + N02 

139  02  + N2 
*140  0,t  + N0  + M 


N02  + NO 

N2  * °2 
0"  + N0 


0 + N2  + M 
0*  + 0o 


02  + W + M 
02  *■  N 
02  ► NO 

o2  > no2 
o2  + no2 

°2  * N2 


02  + N2  + M 


*141 

02(N2)  + M 

-*•  0 

*142 

0*  * 02  + M 

•>  0, 

*143 

Oj  * M 

- o; 

-»•  H30  + NO  + 02 

- N0+(W)2  + M 

-*■  N0*(W)  + W + M* 

- N0+(W)3  + M 

-*•  no+(w)2  + w + m 

- H+(W)3  + HN02 

- N0+(C02)  + n2 
•>  N0+(W)  + N2 

- N0+  + N02 

- °2  + N2 

- N0+  + N 

->  N0+  + N + M 

- o2  + 0 

02(W)  1-  M 

- N0+  0 


NO  + 02 
->  NO*  + 02 

> NO*  + 03 

> NO*  + NO 

- 02(N2)  + M 

-*•  o*  + n2  + f 
-*  oj  + M 

- oj  ♦ o2  + f 


5.0(-10) 
l.l(-27)  -2.0 

5.9(-  2)  -2.0 

1 . 6 ( - 27 ) -2.0 

2 . 2 ( - 2)  -2.0 
7 . 0( - 1 1 ) 

1.0(-  9) 

1.0(-  9) 
2.9(-10) 

5 . 0 ( - 11 ) -0.8 

1 . 2 ( - 1 2 ) -1.0 
6.0(-29)  -2.0 

2 . 0 ( - 1 1 ) -0.4 

2 . 8( -28)  -2.0 
1 . 2 ( - 10 ) 

4 . 5 ( - 10 ) 

6 . 6( - 10) 

0 

1 . 0 ( - 1 6 ) 

0 

0 

3 . 9( -30)  -3.2 

2 . 6 ( - 5)  -3.2 
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Table  A-3.  (Continued). 


Rate  Coefficients 


144 

O^W)  + W 

- h3o+  + oh  + o2 

2 .0(-10) 

145 

O^W)  + W 

- h3o+(oh)  + o2 

1.0(-  9) 

146 

OgCW)  + NO 

-*•  N0+  + W + 02 

l.O(-lO) 

147 

oJdO  + o^1*) 

-*•  o2  + w + o2 

l.O(-lO) 

*148 

oJ(n2)  + 02 

- °;  + n2 

0 

149 

oj  + W 

- oJ(W)  + 02 

1 - 5( - 9) 

150 

o2(w)  + 02 

-oj  * W 

2 . 0( -10) 

151 

0*  + NO 

4 

- NO+  + 02  + 02 

5.0(-10) 

*152 

oj  ♦ N02 

- NO*  + 02  + 02 

0 

153 

o\  * 0 

- 02  + 03 

3.0(-10) 

154 

H+(W)3  + e 

-*■  H + 3(W) 

5.1(-  6) 

155 

H+(W)4  + e 

-*■  H + 4(W) 

6.1(-  6) 

156 

H+(W)5  + e 

- H + 5(W) 

7.4(-  6) 

157 

N0+  + e 

- |n(2D)  + + 0 

4.0(-  7) 

-1.0 

158 

N0+(W)2  + e 

- NO  + 2(W) 

3.0(-  6) 

-0.2 

*159 

N0+(W)3  + e 

- NO  + 3(W) 

5.0(-  6) 

-0.2 

*160 

N0+(N2)  + e 

* NO  + N2 

1 - 5( - 6) 

-1.0 

161 

N2  + e 

- N(2D)  + N 

2 . 7 ( - 7)  . 

162 

02  + e 

-*■  O^A)  + 0 

2. 1(-  7) 

-0.7 

163 

02(W)  + e 

- o2  + w 

1 *5{-  6) 

-0.7 

*164 

oj  + e 

" °2  * °2 

1 • 5 ( - 6) 

-0.7 

165 

e + N02 

> NO" 

4.0(-ll) 
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Table  A-3.  (Continued). 


Reactions 

Rate  Coefficients 

ABC 

166  e + 02  + N2  - 02  + N2 

1 - 0( -31 ) 

167  02  + N2  - e + 02  + N2 

1.9(-12)  +1.5  4990 

168  e + 02  + 02  - 02  + 02 

1.4(-29)  -1.0  600 

169  02  + °2  -*■  e + 02  + 02 

1 . 9(  - 10)  +0.5  5600 

170  e + 03  -*■  0“  + 02 

9 . 0(  - 12 ) +1.5 

171  O'  + N -*■  e + NO 

2.0C-10) 

*172  0"  + NO  -*■  e + N02 

2 . 5( - 10)  -0.8 

173  O'  + 0 -*■  e + 02 

2.0(-10) 

174  0"  + 02( JA)  - e + 03 

3 . 0( - 10) 

175  Oj  + N -►  e + N02 

3.0(-10) 

*176  02  + 0 ■*  e + 03 

1.5(-10) 

177  02  + 02(1A)  ->  e + 02  + 02 

2.0(-10) 

178  0 + e -*■  0"  + hv 

1 . 3( - 15) 

179  e + 02  + H20  ->  02  + H20 

1 . 4 ( -29 ) 

Negative  Ion  Reactions 

*200  CO'  + W + M - COj(W)  + M 

0 

*201  CO'  + NO  N02  + C02 

1 . 8( - 1 1 ) -0.5 

*202  COj  + N02  - NOj  + C02 

2 . 0( - 10) 

203  COj  + 0 - °2  + C02 

l.l(-lO) 

*204  COj(lJ)  + NO  - NOj(W)  + C02 

0 

*205  CO'  + W + M -*•  CO^(W)  + M 

0 

206  CO^  + NO  - NOj  + C02 

4.8(  11) 
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Table  A-3. 


Reactions 


(Continued). 

Rate  Coefficients 


ABC 


207 

CO'  + 0 

- C03  + °2 

1 . 5( - 10) 

208 

c°-  + o3 

O3  + C02  + 02 

1 . 3( -10) 

*209 

COj(W)  + NO 

- NO"(W)  + C02 

0 

*210 

NO"  + W + M 

- NO"(W)  + M 

0 

211 

NO"  + no2 

- NO"  + NO 

0 

212 

NO'  + 03 

- N03  + °2 

1 - 8( — 1 1 ) 

*213 

NO"(W)  + 03 

- NO"(W)  + 02 

0 

*214 

NO'  + W + M 

- NO"(W)  + M 

0 

215 

O’  + C02  + M 

♦ CO3  ♦ M 

3 . 1 ( -28) 

-1.0 

216 

o"  + no2 

- NO"’  + 0 

l-2(-9) 

217 

0"  + 02  + M 

-O'  +M 

l.l(-30) 

-1.0 

218 

0"  + o3 

- °-  + 0 

5 . 3 ( - 10 ) 

219 

O'  + C02  + M 

CO"  + M 

2 . 0( -29) 

-1.0 

*220 

O'  + W + M 

■*  0"(W)  + M 

0 

221 

0'  + N 

-*•  0"  + NO 

l.O(-lO) 

222 

o'  + no2 

- NO"  + 02 

8. 0(-10) 

*223 

O'  ♦ 0 

^ 0"  + 02 

1.5(-10) 

224 

O'  + 02  + M 

°4  + M 

3.5(-31) 

-1.0 

*225 

O'  + M 

- 02  + 02  + M 

2.0(-5) 

-1.0 

226 

°2  + °3 

- °3  + °2 

4 . 0( - 10) 

*227 

O'(W)  + co2 

- CO^  + W 

0 

*228 

0‘(w)  + NO 

* NO"  + H20 

0 

Table  A-3.  (Continued). 


Reactions 

Rate  Coefficients 

ABC 

*229  0"(W)  + 03  - O3  + W + 02 

0 

230  0"  + C02  - CO"  + 02 

5 . 5 ( - 10 ) -0.5 

231  0"  + NO  - NO"  + 02 

l.O(-ll) 

232  O3  + N02  ->  N0‘  + O3 

0 

■•■233  0"  + N02  -*■  NO"  + 02 

2.8(-10) 

*234  O3  + 0 °2  + °2 

2.5(-10) 

235  0^  + C02  - CO"  + 02 

4.3(-10) 

*236  CO"  + 02  ^ 0"  + C02 

8 . 4( - 12 ) 3000 

*237  0“  + W - 0"(W)  + 02 

0 

238  0"  + NO  - NO"  + 02 

2 . 5 ( - 10) 

239  0"  + 0 - 0"  + 02 

4 . 0( - 10) 

*240  0“  + O3  -*•  O3  + 02  + 02 

3.0(-10) 

241  NO"  + NO  -*■  NO"  + N02 

1.5(— 11) 

242  M+  + M"  -*■  Products 

6.0(-8) 

243  M+  + M"  + M •>  Products 

5 . 6( -26)  -1.5 

APPENDIX  B 

IONOSPHERIC  SOURCE  FUNCTIONS 


For  sensitivity  analyses,  we  require  natural  ionization  source 
functior  models  that  incorporate  the  various  physical  relationships 
explicitly.  We  therefore  want  to  avoid  models  for  effective  source 
functions  that  are  derived  from  observed  ionization  or  propagation 
conditions.  W'e  have  taken  models  based  on  those  in  Reference  B-l, 
since  they  are  computationally  simple  and  incorporate  most  of  the  de- 
sired physical  relationships.  We  have  modified  the  functional  form 
in  some  cases  to  make  explicit  certain  relationships. 

GALACTIC  COSMIC  RAY  IONIZATION 

The  ion  production  rate  from  galactic  cosmic  rays  is  computed  from 
qc  = 4(30+  250 f sin40)  * 10_2°(M)  cm"3s_1  (B-l) 

where 

0 = geomagnetic  latitude 

(M)  = total  neutral  particle  concentration  at  deposition 
point  (cm  3) 

f = 4 - 3S  (B-2) 

S - relative  solar  activity  (0  to  1). 

One  measure  of  solar  activity  is  the  sunspot  number  SSN  , which  varies 
from  about  10  for  minimum  conditions  to  about  110  for  maximum  conditions. 
So  the  relative  parameter  S may  be  defined 

S = (SCM  - 10)/ 100  . (B-3) 

oiN 

Reference- B-l . Knapp,  W.S.,  ct  al.  Reaction  Rate,  Collision  Frequency 
and  Ambient  Ionosphere  Models  for  Use  in  Studies  of  Radio  Propaga- 
tion in  a Uuclear  Environment,  66TMP-18  (DASA  1765),  General  Ulcc- 
tric  Company — TliMI’O , March  1966. 


» 
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Note  that  the  cosmic  ray  ion  production  rate  varies  about  a factor 
of  4 with  solar  activity,  and  that  it  varies  in  an  inverse  manner 
with  maximum  ion  production  during  minimum  solar  activity.  The  strong 
latitude  dependence  results  from  deflection  of  the  charged  particles 
by  the  earth's  magnetic  field.  The  given  expression  is  applicable 
only  to  altitudes  above  40  km,  but  as  will  be  shown  later,  there  is 
negligible  IIP  signal  absorption  below  this  altitude. 

SOLAR  LYMAN-ALPHA  IONIZATION 

The  ionization  of  nitric  oxide  (NO)  due  to  solar  Lyman-a  flux  is 
given  by 

qLa(N0+)  = 6,0  * 10'7(NO)fLcexp[-1.0  * 10_2°I02]  (B-4) 


where 


(NO)  = nitric  oxide  concentration  at  deposition  point  (cm""*) 


(B-S) 

IQ2  = line  integral  of  (02)  along  the  path  of  the  solar 
flux  (cm"2) . 


If  we  take  21  percent  of  the  air  mass  as  (02) , and  employ  hydrostatic 
equilibrium,  the  line  integral  may  be  computed  from 


IQ2  = 4.7  x 1018P  secy  cm"2  (B-6) 

_2 

P = atmospheric  pressure  at  deposition  point  (dyne  cm  ) 

X = solar  zenith  angle. 


The  form  secx  is  a flat  earth  approximation  and  will  overestimate 
(underest imatc  q ) for  solar  zenith  angles  greater  than  about  80  de- 
grees. The  relative  solar  activity  S is  from  liquation  B-3. 

SOLAR  X-RAY  IONIZATION 

Solar  x-rays  include  that  part  of  the  solar  spectrum  in  the  0-2QA 
wavelength  band,  in  which  the  energy  transport  is  typically  described 
in  terms  of  mass  absorption  coefficient.  The  ion  production  rate  is 
given  by 
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♦ £0*' 


<*X  = 6p/  FS(X)VX)  exp[-yro(X)m]  dX 


APPENDIX  B 


(B-7) 


where 

B = number  of  ion  pairs  per  unit  absorbed  energy 
(1.8  x 1010  erg-1) 

_3 

p = air  mass  density  at  deposition  point  (g  cm  ) 

X = photon  wavelength  (A) 

F_(X)  = solar  spectrum:  flux  at  top  of  atmosphere 

, -2  -lrl. 

(erg  cm  s A ) 

2 -1 

U (X)  = mass  absorption  coefficient  for  air  (cm  g ) 

ra  -2 
m = air  mass  penetrated  to  deposition  point  (g  cm  ) . 

The  integral  has  been  computed  in  Reference  B-l  using  analytic  expres- 
sions to  fit  the  solar  spectrum  and  mass  absorption  coefficients.  The 
integral  is  performed  piecewise  due  to  natural  breaks  in  the  analytic 
curve  fits  and  also  because  the  solar  intensity  variation  over  the  solar 
cycle  is  a function  of  wavelength.  For  X = 0-  8A  , the  solar  intensity 
is  assumed  to  vary  by  a factor  of  1000;  for  X = 8 - 20 A,  the  variation 
is  about  a factor  of  45.  The  results  of  integration  produce  a function 
of  air  mass  penetrated  (m) : 


where 


qx  = pflxtll(ra)  + I2(m)]  + ^ZxV"0 


f,  = 0.001-  + 0.999  s 

lx 


f.  = 0.022  + 0.978  S 

2x 


(B-8) 


(B-9) 
(B- 10) 


and  S is  defined  in  Equation  B-3.  If  we  take  hydrostatic  equilibrium: 


m = — sccy 
g 


(B-ll) 


and  the  absorption  integrals  may  be  expressed  as  functions  of  P secX. 
These  arc  shown  in  Figure  B-l,  and  for  computational  purposes  arc  fit 
by  the  following  expressions: 
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IA (x)  = 4.80  x 10 


= 6.04  x io°  cxp(-0.23x) 


= 2.71  x io°  cxp(-0.15x) 


= 3.89  x I0ll/x3*46 


I2tx)  = 1.56  x io* 


1 0 T 

= 4.47  x 10  /(2.871  ♦ xJ) 


= 1.20  x lO^/x6-43 


I3(x)  = 1.00  x io‘ 


= 1.09  x io*  exp (-8. 70x) 


= 5.40  x io  exp(-2.60x)/x 


= 2.40  x io1*  exp (-1 . 80x) 


X £ 1 


1 < x s;  io 


( B- 12) 


10  < x s 30 


30  < x 


x £ 0.1 


0.1  < x i 10 


(B-13) 


10  < x 


0.001  < x £ 0.1 


0.1  < x <:  1 


(B-14) 


1 < x 


where 


x = P secx 


SOLAR  SOFT  X-RAY  AND  UV  IONIZATION 


The  first-order  ionization  of  (0  ) , (0  ) , and  (N2)  due  to  soft 


x-rays  and  ultraviolet  energy  is  computed  from 


quv(°+)  = (0)yV0(A)  FS(X)  cxp[-T02(A) -TN2(A)}dX 
‘IUV(023  = (02)/CT02(X)  FS(A)  cxp[-T0,(X)  " TN2(A)]d^ 
Muv(n^  = (n2)/"°n2(A)  FS(A)  CXP[*T02(A)  ‘TN2(A)ldA 


(B-1S) 


(B- 16) 


(B-17) 


where 


(0)  = atomic  oxygen  concentration  (cm’  ) 


(O  ) = molecular  oxygen  concentration  (cm  ) 


(N2)  = molecular  nitrogen  concentration  (cm  ) 





r-  +*-**"- 


! 

i 


°0(X) 

°02^ 

°i«<X> 

FS(X) 


X02^ 

TN2<X) 

°02(X) 


2 

ionization  cross  section  for  0 (era  ) 

2 

ionization  cross  section  for  (cm  ) 

2 

ionization  cross  section  for  ^ (cm  ) 

solar  photon  flux  spectrum  at  top  of  atmosphere 
(photons  cm  2 . 

°02W  f”° 2)dl 

°N2(‘,/”N2)dL 

•'z 

2 

total  absorption  cross  section  for  0^  (cm  ) 

2 

total  absorption  cross  section  for  Nj  (cm  ) 


( B- 18) 
(B- 19) 


The  line  integrals  defining  the  optical  depths  Tq2  and  (Equations 
B-18  and  B-19)  are  along  the  path  of  solar  flux  down  to  altitude  z. 

The  ionization  and  absorption  cross  sections  arc  quite  structured,  and 
since  the  photon  flux  spectrum  is  not  sufficiently  well  established 
or  constant,  broad-band  averages  arc  employed  for  both  the  flux  and 
cross-sections.  These  are  listed  in  Table  B-l.  Ke  compute  the  line 
integrals  of  liquations  B-18  and  B-19  from 


f f<pP 

J (°2)dL  = 52^  SCCX 

z 

rr ...  fN2p 

J °VdL  = 28 SCCX 


(B-20) 

(B-21) 


where  f„  and  f,,_  arc  the  local  volume  mixing  ratios  for  (0.)  and  (N_) , 
respectively,  and  tn^  = 1.66  * 10"  g . So  that  the  integrals  of  F.qua- 
tions  B- 15  through  B-17  may  be  precomputed  as  a function  of  a single 
variable,  we  make  the  following  approximations  for  the  volume  mixing 
ratios 


78 


APPENDIX 


Table  0-1.  Ionization  and  absorption  cross  sections  (units  of  10"^ 
cm2)  and  free  space  solar  fluxes  from  10  to  1026 A . 


A (A) 

1 0”**  <P0A 

0(0) 

o'  <Oa) 

o'  (N2) 

0(0,) 

o (Na) 

1025.7 

26 

1.0 

1.7 

1000-1027 

15 

1.3 

1.8 

989.8 

5 

1.8 

2.2 

0.4 

977.0 

30 

3 

3.7 

0.8 

972.5 

10 

25 

30 

280 

949.7 

5 

5 

5.6 

1.9 

911-1000 

37 

9 

11 

10 

850-911 

95 

3 

9 

11 

12 

796-850 

25 

3.5 

18 

22 

10 

700-796 

50 

5 

25 

21 

30 

26 

600-700 

47 

10 

25 

26 

30 

33 

584.3 

29 

13 

25 

24 

30 

30 

500-600 

30 

15 

25 

20 

30 

25 

400-500 

24 

11 

25 

20 

30 

25 

303.8 

43 

10 

15 

9 

17 

10 

300-400 

29 

10 

15 

12 

16 

15 

230-300 

31 

9 

13 

8 

15 

9 

170-230 

33 

7 

10 

5 

11 

6 

110-170 

3.5 

2.8 

5 

1.8 

5 

1.8 

80-110 

2.4 

?.! 

2 

0.5 

2 

0.5 

60-80 

1.8 

0.5 

1 

0.25 

l 

0.25 

30-60 

1.5 

0.25 

0.5 

0.12 

0.5 

0.12 

20-30 

0.12 

0.06 

0.1 

1.0 

0.1 

1.0 

10-20 

0.02 

0.30 

0.6 

0.3 

0.6 

0.3 

.1 


f02  - °'21 


P £ 1.65 


= 0.21  - 0.017  in  (1.65/P)  , P < 1.65 


(B-22) 


fN2  = °-79 


P £ 1.65 


= 0.79  - 0.012  in  (1.65/P)  , P < 1.65 


(B-23) 


The  ion  production  rates  arc  then  given  by 


quv(0  > = pf0fuvI4(p  sinx) 


(B-24) 


Muv(°  ) = pfyv1^’  sinx) 


(B-25) 


clUV(N2)  = pfUVr6(P  SlnX) 


(B-26) 


where  fQ  is  the  volume  mixing  ratio  to  atomic  oxygen  and 


fuv  = i(1  + s)  * 


(B— 27) 


and  the  relative  solar  activity  S is  from  Equation  B-3.  The  functions 

I.,  I_,  and  I,  arc  illustrated  in  Figure  B-2.  For  computational  pur- 
4 D o 

poses,  we  employ  the  following  curve-fit  expressions: 


I4 (X)  = 2.94  x io14  exp(-129X)/X°-432 
= 3.20  x IQ7  exp(+28.7X)/X3-63 


x £ 0.0316 


0.0316  < x £ 0.1 


= 9.23  x 10llexp(-7.47X)/X°'741 


I5(X)  = 1.24  x io1J  exp(-198X) 

= 6.16  x io10  cxp(+14.6X)/X1-69 


0.1  < x (B-28) 


x £ 0.01 


0.01  < x s 0.1 


j 


= 2.63  x lo1J  cxp(-7.04x) 


, 0.1  < x . (B-29) 


0.0316  < x * 0.1 


0.1  < x 


( B - 30) 


I6(X)  = 5.70  x 1013  cxp(-126X)/X0,667  , x i 0.0316 

=■  1.36  x io7  exp(+25.9X)/X3'69  , 0.0316  < x*  0.1 

= 3. 54  x io11  exp(-7.45X)/X0,723  , 0.1  < x (B-30) 

FINAL  EXPRESSIONS  FOR  ION  PRODUCTION  RATE 

The  particular  ions  produced  by  the  UV  source  are  0 , 02,  and  N2> 
as  noted  in  liquations  B-24  through  B-26*.  The  Lyinan-alpha  source  (liqua- 
tion B-4)  produces  the  NO*  ion.  Most  of  the  ionization  produced  by 
the  cosmic  ray  and  X-ray  sources  results  from  the  cascade  process  of 
secondary  and  tertiary  energetic  electrons.  For  these  sources,  the 
following  partitioning  of  the  ionization  is  employed,  based  on  calcu- 
lations in  Reference  B-2: 


q(N+)  = (qc  + qx)[UN(N)  *0.20nN2(N2)  /T 
q(0+)  = (qc*qxuyo)  + 0.27n02(O2)]/T 
q(N*)  = (qc+qx)[0.8nN2(N2))/T 
q(02)  = (qc+qx)[0.73n02(02)]/T 


(B— 31) 


(B-32) 


(B-33) 


(B-34) 


q (NO  )=  (qc  + qx)[nNO(NO)]/T 


(B-35) 


where 


nN  = 0.54 
n0  = 0.67 

°N2  = 1-°° 
n02  ^ 
\'0=  l*45 


t = nN(N)  ♦ n0(0)  ♦ nN2(N2)  ♦ nQ2(02)  ♦ nNQ(NO)  . 

Reference  B-2.  Stephens,  T.L.,  and  A. L.  Klein,  Electron  Energy  Deposi- 
tion in  the  Atmosphere , GU75TMP-7,  General  Electric  Company — TEMPO, 
May  1975. 
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1 


Note  that  for  typical  mesospheric  mixing  ratios  (78  percent  N2  > 21 
percent  0 ^ » less  than  0.1  percent  0,  less  than  one  part  per  million 
NO),  the  ion  partition  of  the  cosmic  and  X-ray  sources  is  as  follows: 


ION 

FRACTION 

N* 

.154 

0* 

.063 

.624 

.169 

N0+  (less  than  about  10  6) . 

ttu.s.  GOVERNMENT  PRINTING  OFFICE.  1979  677-01 7A« 
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